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Abstract — We consider the performance modeling and evalua- 
tion of network systems powered with renewable energy sources 
such as solar and wind energy. Such energy sources largely 
depend on environmental conditions, which are hard to predict 
accurately. As such, it may only make sense to require the 
network systems to support a soft quality of service (QoS) 
guarantee, i.e., to guarantee a service requirement with a certain 
high probability. In this paper, we intend to build a solid 
mathematical foundation to help better understand the stochastic 
energy constraint and the inherent correlation between QoS and 
the uncertain energy supply. We utilize a calculus approach 
to model the cumulative amount of charged energy and the 
cumulative amount of consumed energy. We derive upper and 
lower bounds on the remaining energy level based on a stochastic 
energy charging rate and a stochastic energy discharging rate. 
By building the bridge between energy consumption and task 
execution (i.e., service), we study the QoS guarantee under the 
constraint of uncertain energy sources. We further show how 
performance bounds can be improved if some strong assumptions 
can be made. 

Index Terms — Stochastic Network Calculus, Performance Eval- 
uation, Renewable Energy, Energy Scheduling 

I. Introduction 

In the last decades, there have been increasing demands on 
computing, communication and storage systems. The strong 
demands drive modern IT infrastructures to extend and scale 
at an unprecedented speed, raising serious "green" related con- 
cerns about high energy consumption and greenhouse emis- 
sion. Addressing this problem, the use of renewable energy 
sources, such as solar and wind energy, plays an important role 
in the support of sustainable computing [5|. Various energy 
harvesting devices have been developed and broadly used 
in real-world applications. For example, to support perpetual 
environmental monitoring, solar energy has been used to 
power tiny sensor nodes deployed in the wilderness JS). 

While the benefit of using renewable energy is clear, its 
application can pose great challenges in terms of assuring a 
satisfactory quality of services (QoS) for computing systems. 
First, renewable energy sources are generally unreliable and 
hard to predict. The uncertainty in the energy supply makes 
it extremely hard to support strict QoS requirements. As a 
result, it only makes sense to demand the system to support 
a soft QoS guarantee, i.e., to guarantee the QoS requirement 
(e.g., delay, throughput) with a certain probability. Second, 
although modern computing and communication devices could 
be equipped with rate-adaptive capabilities ||23ll , allowing 
the devices to "smartly" schedule the execution of various 
tasks, it is nevertheless non-trivial to design good scheduling 



algorithms that can effectively utilize a limited and variable 
energy resource. Many aspects come into play, including for 
example, the predicted energy charging rate, the tolerable 
range of a QoS guarantee, the rate-adaptive features of the 
device, and so on. Third, QoS support and task scheduling 
become even harder when the network includes multiple nodes 
powered by renewable energy, because the task scheduling and 
energy management of one node may have direct impact on the 
other nodes. Such dependency makes QoS support extremely 
hard. 

Given these challenges, we perceive the strong need for 
a generic analytical framework for performance modeling 
and evaluation of a network system using renewable energy 
sources. Such a theoretical model should (1) capture the 
stochastic features in energy replenishment and energy con- 
sumption, (2) provide good guidelines on the schedulability 
of given tasks under uncertain energy constraints, (3) be 
applicable to a large group of systems where the power-rate 
function fTTll . l|23l may be of different forms, and (4) be able 
to analyze a single node as well as a network system. 

Various methods have been developed for task schedul- 
ing and performance analysis of energy constrained system, 
including for example Markov chain based methods [16|, 
[20 1, calculus approaches ITD . Il23l . prediction-based ap- 
proaches |7|. [21], and so on. Nevertheless, no analytical 
model so far is sufficient to meet all the above requirements. 
We are thus motivated to develop a more generic analytical 
framework to fill the vacancy. In this paper, we make the 
following contributions: 

> We build a theoretical model for performance evaluation 
of network systems with renewable energy sources. Our 
model is based on recent progress in stochastic network 
calculus, but it much extends the concepts of the tra- 
ditional stochastic network calculus by introducing en- 
ergy charging/discharging models. Such extension is non- 
trivial since the new concepts require special treatment in 
the derivation of performance bounds. 

• We derive the stochastic upper and lower bounds on 
a node's residual energy level. These bounds provide 
fundamental guidance in the energy management of a 
system with renewable energy supply. 

• By using a generic power-rate function to bridge the 
task execution and its energy consumption, we derive 
the stochastic performance bounds on delay and system 
backlog. We study both the single-node case and the 
network case. 
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• We extend the model to analyze systems that utilize 
multiple energy sources. 

• We point out a method to improve the performance 
bounds if some strong assumptions, such as independence 
between multiple energy sources, can be made. 

The rest of the paper is organized as follows. We introduce 
the basic notations in SectionUI] In SectionUm we present new 
energy models to capture the stochastic features in the energy 
charging and discharging processes. Based on the new energy 
models, we derive the stochastic bounds on a node's remaining 
energy. We analyze the performance bounds of a single node 
with respect to delay and system backlog in Section [IV] The 
performance of a network system is provided in Section [V] 
Methods of handling multiple energy sources and improving 
performance bounds are introduced in Sections [VT] and IVIII 
respectively. Related work is discussed in Section IVIIII The 
paper is concluded in Section [IX] 

II. Basic Notations 

We first introduce the basic notations following the conven- 
tion of stochastic network calculus (9), 1101 . 1131 . We denote 
by T the set of non-negative, wide-sense increasing functions, 
i.e., 

T = {/(•) : VO < x < y, < f(x) < f(y)}, 

and by T the set of non-negative, wide-sense decreasing 
functions, i.e., 

T = {/(•) : VO < x < y, < f(y) < f(x)}. 

For any random variable X, its distribution function, denoted 
by 

F x {x) = Prob{X < x}, 

belongs to J 7 , and its complementary distribution function, 
denoted by 

F x (x) = Prob{X > x}, 

belongs to T. 

For any function f(t), we use f'(t) to denote its derivative, 
if it exists. 

The following operations will be used in this paper: 

• The (min, +) convolution of functions / and g under the 
(min, +) algebra HI, j4), lfl2ll is defined as: 

tf<8>g)(t)= taf ms) + g(t-s)}, (1) 

0<s<t 

• The (min, +) deconvolution of functions / and g is 
defined as: 

(f<2>g)(t) = su V {f(t+s)-g(s)}. (2) 

• The (max, +) convolution of functions / and g is defined 

as: 

(m)(t)= sup {f( S ) + g(t-s)}. (3) 

0<s<t 

• The pointwise infimum or pointwise minimum of functions 
/ and g is defined as 

(fAg)(t)=min[f(t),g(t)}. (4) 

In addition, we adopt: 



• [x] + = max{x, 0}, 
« [x]i = min{x, 1}. 

Throughout this paper, we assume that the energy charging 
curve and the energy discharging curve, both of which will 
be defined later, are non-negative and wide-sense increasing 
functions. In this paper, C(t) and C*(t) are used to denote the 
cumulative energy amount that has been charged and depleted 
in the time interval (0, t], respectively. For any < s < t, 
let C{s,t) = C(t) - C(s) and C*(s,t) = C*(t) - C*(s). 
Similarly, A(t) and A*(t) are used to denote the cumulative 
amount of data traffic that has arrived and departed in time 
interval (0,i], respectively, and S(t) is used to denote the 
cumulative amount of service provided by the system in time 
interval (0,t], For any < s < t, let A(s,t) = A(t) - 
A(s),A*(s,t) = A* (t)- A* {s), and S(s,t) = S(t)-S(s). By 
default, A(0) = A*(0) = 5(0) = 0, and C(0) = C*(0) = 0. 

III. Stochastic Energy Models and Residual 
Energy 

A. Stochastic Energy Charging and Discharging Models 

Definition 1: The stochastic energy charging model: The 

cumulative energy amount C(t) is said to follow stochastic 
energy charging (s.e.c.) curves ol\ £ T and 012 € T with 
bounding functions fx £ T and /2 G T, respectively, denoted 
by 

C ~sec< fl, Oil, hi &2 >, 

if for all t > s > and all x > 0, 

Prob{ inf [C(s,t)-a 1 (t-s)>x]}>f 1 (x), (5) 

0<s<t 

and 

Prob{ sup [C(s,i) ~a 2 (t- s) > x]} < f 2 (x). (6) 

0<s<t 

We call the curves ct\ and a 2 the lower curve and the 
upper curve, respectively, and the functions /1 and h me 
lower bounding function and the upper bounding function, 
respectively. 

Remark 1: The practical meaning of (0 is to lower bound 
the cumulative amount of charged energy, i.e., the energy 
harvester should provide a minimal level of energy with a high 
probability. The meaning of © is that the cumulative amount 
of charged energy may be upper bounded. For instance, in the 
case of solar-powered sensor nodes, the total charged energy 
should not beyond the best case scenario, e.g., sunny all the 
time. 

Remark 2: The (p, cr\, <72)-source in the harvesting theory 
ifm is a special case of our stochastic energy charging model, 
by defining ai(t) — pt — a%, a 2 (t) = pt + a 2 , and setting the 
bounding functions to 0. 

Definition 2: The stochastic energy discharging model: 
Under the constraint of charged energy amount C(t), the 
cumulative discharged energy amount C*(t) is said to have 
stochastic energy discharging (s.e.d.) curves ft £ J and 
j3 2 G T with bounding functions g\ G T and g 2 G T, 
respectively, denoted by 

C* ~sed< gi,Pi,g 2 ,& 2 >, 
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if for all t > and all x > 0, 

Prob{C®p 1 (t)-C*(t)>x}>gi(x), (7) 

and 

Prob{C gi /3 a (t) - C* (t) > x} < ff2 (x) . (8) 

We call curves j3\ and 02 the upper curve and the lower curve, 
respectively, and the functions 51,32 the upper bounding 
function and the lower bounding function, respectively. 

Remark 3: The practical meaning of the above definition is 
as follows. If we use /3i(t) to denote the cumulative virtual 
energy discharging amount (i.e., the energy budget that a 
service should follow), and if we want to use (3\(t) to upper 
bound the actual discharged energy C*(t), then under the 
constraint of cumulative charged energy amount C(t), at any 
time instance t, we have 

C*{t)<^{s) + C{t-s). (9) 

Note that we should not replace /3i(s) with C*(s) since 
otherwise we obtain C*(t — s) < C(t — s) for any time 
interval (s,t], which is a constraint too restrictive. Because 
the inequality (0 holds for any s < t, we have C*(t) < 
inf < s < t {/3 1 (s)+C(t-s)}, which is C*(t) < C®Pi{t) by the 
definition of (min,+) convolution. Inequality (0 in Definition|2] 
represents the stochastic version of the above relationship. In 
practice, we may also want to lower bound the discharged 
energy (with a curve f5 2 ) to effectively utilize the harvested 
energy, since overcharged energy beyond the energy storage 
capacity will be wasted. In this case, the stochastic lower 
bound is represented by the inequality ©. 

Since C(t) and C*(t) denote the cumulative amount of 
charged energy and the cumulative amount of discharged 
energy, respectively, the remaining energy amount, denoted 
by E{t), in the system at time t can be calculated as: 

E(t) = C(t)-C*(t). (10) 

B. Stochastic Analysis on Remaining Energy 

For any energy-aware scheduling, we first need to under- 
stand the feasible energy depletion region. Particularly, we 
need to answer the following question: what are the lower 
and upper bounds of Prob{E{t) > x}? 

The practical meaning of the lower bound comes from 
reliability consideration: the remaining energy at any time 
instance should be with a high probability larger than a 
minimal threshold value (called safety threshold) to guarantee 
the reliable operation of the system. The upper bound is related 
to the effective use of the energy source: since the capacity of 
energy storage is limited (in this case, x could be considered 
as the upper limit of the energy capacity), tasks should be 
scheduled to consume energy effectively so that the probability 
that the system is overcharged is small. 

To answer the above question, we start with the following 
lemma. 

Lemma 1: Assume that the complementary cumulative dis- 
tribution functions ( CCDF) of random variables X, Y are 
Fx £ J~ and Fy 6 T, respectively. Assume that fi{x) < 
Fx(x) < f 2 {x) and g\(x) < Fy(x) < 32 (%)■ Denote 



Z = X+Y. Then no matter whether X and Y are independent 
or not, there holds for Vx > 0, 

F z (x) < f 2 ®92(x), (11) 

and 

F z {x) > fi®gi{x) - 1. (12) 

The proofs of all lemmas are referred to Appendix of this 
paper. 

Theorem 1: Assume that the system has an energy charging 
source with cumulative charged energy amount C(t) ^ sec < 
fx, Q!i, /2, 012 > and it provides service with cumulative 
depleted energy amount C*(t) ~ se d< gi,Pu92,fh >■ The 
remaining energy amount of the system at any time instant t, 
E(t), is lower bounded by 

Prob{E(t) > x}> h®gi{x -a x ft(0)) - 1, (13) 

and is upper bounded by 

Prob{E{t) > x} < / a ® g 2 (x - a 2 /3 2 (0)), (14) 

Proof: We only prove the lower bound, since the upper 
bound can be easily proved following the same argument using 
the curves a 2 and f3 2 , the bounding functions f 2 and g 2 , and 
Inequality ( fTTT i. 

For any t > s > 0, we have 

E(t) =C(t)-C*(t) 

=C(t) -C®px(t) + C® ft (i) - C*(t) 

= sup {C(M) - cti(t - s) + ai{t - s) - ft (t - s)} 

0<s<t 

+ C®l3 1 (t)-C*(t) 
> sup { inf {C(s,t) - at(t - s)} 

0<s<t 0<s<t 

+ a 1 (t-s)-pi(t-s)} + C®[3 1 (t)-C*(t) 
= inf {C(s, t) — ai(t — s)} + sup {aii(t — s) 

0<s<t 0<s<t 

-j8i(t-«)} + C7®|8i(i)-C*(t) 
= inf {C(s,t)-ai(t-s)} + sup{a 1 (t)-/3 1 (t)} 

0<s<t t>0 

+ C®px{t)-C*{t) 
= inf {C{s,t)-a x {t- s)} + C®p x {t)-C*{t) 

0<s<t 

+ ai 0i(O) 

(15) 

If we consider mf < s <t{C(s, t) — ai(t — s)} and C®j3i{t) — 
C* (t) as two random variables depending on t, Inequality (fT3l l 
follows according to the definitions of s.e.c. and s.e.d. curves 
and Inequality (fT2l in Lemma [TJ ■ 
Remark 4: Theorem [TJ provides the following fundamental 
guidance in the energy management of systems powered with 
renewable energy: 

• Given the statistical feature of a renewable energy source 
and a safety threshold on remaining energy level, we 
can find a suitable energy consumption budget to service 
requirements (i.e., /3i(i) in Theorem [TJ based on dT3b , 
so that the system has sufficient energy with a certain 
probability. 
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• Given the statistical feature of a renewable energy source 
and the maximum capacity of energy storage in the 
system, we can obtain a minimal energy expenditure rate 
based on (fl4l >. such that the probability that the system 
is overcharged is lower than a certain probability. 

• Overall, given a renewable energy source, a safety thresh- 
old energy level, and an upper limit for energy storage, 
the two inequalities in Theorem Q] describe the stochasti- 
cally feasible region of energy budget for services. 

IV. Stochastic Service Guarantee of A Single 
Node 

A. System Model 

With the above energy models, we next study the problem 
of providing a stochastic service guarantee given an uncertain 
energy supply. We start with the single-node case. 

For the paper to be self-contained, we need to briefly 
introduce the core concepts in traditional stochastic net- 
work calculus- stochastic traffic arrival curves and service 
curves |10|. For traffic arrivals, we have the following model. 

Definition 3: The v.b.c. model: A flow A(t) is said to 
have a virtual-backlog-centric (v.b.c.) stochastic arrival curve 
a G J- with bounding function / € T, denoted by 

A ~„&< /, a >, 

if for all t > and all x > 0, it holds (9), iflOl 

Prob{ sup [A(s, t) - a(t - s)] > x} < f(x). (16) 

0<s<i 

For service models, we have the followings. 

Definition 4: The s.c. model: A server is said to provide 
a flow A(t) with a stochastic service curve (s.c.) (3 G T with 
bounding function g G T, denoted by 

S ~ sc < g,P >, 

if for all f > and all x > 0, it holds J9), IflOl 

Prob{ sup [A ® /3(s) - A*(s)} > x} < g{x). (17) 

0<s<t 

Remark 5: The s.c. model is adopted in this paper for 
ease of expressing the results, particularly the concatenation 
property for the network case analysis in Section [V] However, 
the s.c. model may be too restrictive. In the literature, a 
variation of the model is available, which is called weak 
stochastic service curve (w.s) in The w.s. model is much 
less restrictive and has been widely adopted in the stochastic 
network calculus literature. We would like to stress that the 
w.s. model can also be used here, and it can be verified (e.g. 
see 1 9 1) that the delay and backlog bound results in this paper 
remain unchanged with the w.s. model. The affected one is 
the analysis on the network case, which would have much 
complicated expression. 

Definition 5: The s.s.c. model: A server is said to provide 
a flow A(t) with a strict stochastic service curve (s.s.c.) [3 G T 
with bounding function g G T, denoted by 

S ~ ssc < g,0 >, 



if during any period (s, t] and for any x > 0, it holds (9), IflOl 

Prob{S(s, t) < /3(t -s)-x}< g(x). (18) 

Remark 6: The s.s.c. model is to decouple the service and 
traffic arrivals. It is useful in energy outage analysis as shown 
later in Section IIV-CI 

The following measures are of interest in service guarantee 
analysis under network calculus: 

• The backlog B(t) in the system at time t is defined as: 

B(t) = A(t) - A*(t). (19) 
> The delay D(t) at time t is defined as: 

D(t) = inf{r > : A(t) < A*(t + r)}. (20) 

B. Stochastic Service Guarantee 

To ease description, we define the following terms: 
Definition 6: Power-rate function 11231 is a function that 
translates the amount of service to the amount of consumed 
energy. We use a generic notation V to denote the power-rate 
function. 

Definition 7: Energy-oblivious service curve of a system 
is the service curve that the system would have if there were 
no energy constraint. 

Note that most modern devices, such as computer servers 
and wireless transceivers, have the capability of adjusting 
processing/transmission rate. Associated with a rate, there 
is a corresponding power expenditure that is governed by 
the power rate function. Generally speaking, a low process- 
ing/transmission rate requires low energy consumption. The 
power-rate function is system dependent. For example, for 
most encoding schemes in wireless communication, the re- 
quired power is a convex function of the rate 11231 . In addition, 
there are some devices, e.g., Atmel microprocessors, whose 
energy consumption is dominated by its on state. That is, as 
long as the device is powered on, it consumes energy at a 
roughly constant speed, no matter whether or not its CPU 
remains idle or executes tasks. To avoid those different details, 
we use a generic notation V to build a broadly-applicable 
model. 

Given a cumulative energy amount C(t), the service amount 
that C(t) can support, denoted as Sc(t), can thus be calculated 

as: t 

Sc(t) = I p-^C'ix^dx, (21) 
Jo 

Similarly, given a service amount S(t), the corresponding 
energy consumption amount, denoted by Cg, can be calculated 

as: t 

C* s (t) = [ V(S'(x))dx. (22) 
Jo 

Remark 7: We assume that the inverse of power rate func- 
tion V, the derivative of the energy charging curve C, and the 
derivative of the service curve S, all exist. This assumption has 
been used and justified in [23 1. This assumption is reasonable, 
because (1) there is usually a one-to-one mapping between 
the served data amount and the used energy amount, and (2) 
cumulative energy charging amount and cumulative service 
amount are mostly continuous. Nevertheless, even if the above 
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assumption may not be true in some specific situations, 
numerical approximation can be used to estimate Sc and Cg. 
In the rest of the paper, we will use Equations ( f2TT > and (1221 
without giving further explanation. 

With all notations being introduced, we are ready to answer 
the following question regarding a QoS guarantee: 

Assume that a data flow A(t) with the traffic arrival 
curve A ~„&< /, a > is input into a node, which has 
an energy-oblivious service curve S ^ S c< g,P >■ Assume 
that the system is powered by an energy source following 
C{t) r ^ S ec< fi, oti, fz> ot 2 >■ For an y time t, what are the 
stochastic bounds on the flow's delay D(t) and backlog B(t)? 

The main difficulty in solving the above question is that the 
data flow may not be able to get a service following the service 
curve f3 due to the energy constraint. In other words, a service 
is possible only if the system has enough energy. By "enough 
energy", we mean that the system's energy level is above 
the safety threshold. To simplify presentation, we assume the 
safety threshold value is in this section. Otherwise, trivial 
modification is required for the following performance results. 

We have the following important technical lemmas. 

Lemma 2: Assume that functions /i,/2,/3 G T, and as- 
sume that /i(0) = / 2 (0) = / 3 (0) = 0. For any t > 0, it holds 
that 

fx ® (h A f 3 )(t) < (/i ® f 2 (t)) A f 3 (t) (23) 
Lemma 3: If Xi > X 2 > and X 3 > X 4 > 0, it holds 



Xi A X 4 - X 2 A X 3 < Xi - X 2 + X 3 - X 4 . (24) 

Lemma 4: | 1 1 Consider a random variable X. For any x > 
0, Prob{[X]+ >x} = Prob{X > x}. 

Theorem 2: Assume that a node has an energy-oblivious 
service curve S ^ S c< g, P >■ Assume that the node 
is powered with an energy source following C(t) ~ sec < 
/li^ij fi-,012 >■ The actual service available to the task, 
denoted by S e (t), follows S e (t) ^ sc < g,/3 A d 2 >> where 



g(x) = g®f 2 (x), 

d 2 (t) = / V- 1 (a 2 (x))dx. 
Jo 

Proof: It is clear that at any time t > 0, the actual 
provided by the system. 



S e (t) 



S(t) if S(t) < S c (t) 
Sc (t) otherwise, 



(25) 
(26) 
service 

(27) 



Note that in the above, the first inequality is due to 
LemmalU the second equality is because that A*(s) = S e (s), 
and the last inequality is based on Lemma [3] 

We thus have 

sup {A (g> (/3 A d 2 )(s) - A*(s)} 

0<s<t 

< sup {[A ® /?(s) - A*(s)}+ + [S c {s) - d 2 (s)] + } 

0<s<t 

< sup {[A®p(s)-A*(s)}+}+ sup {[S c (s) ~d 2 (s)} + } 

0<s<t 0<s<t 

(29) 

Based on the definitions of the s.c. curve and the s.e.c. curve, 
and Lemma |U we have for any x > 0, < s < t, 

Prob{ sup [A ® I3(s) - A*{s)] + > x} < g(x), (30) 

0<s<t 

Prob{ sup [S c {s) - ci 2 (s)} + >x}< f 2 (x). (31) 

0<s<t 

For OTT i. we remind the readers that the calculations of Sc 
and d 2 do not change the bounding function. The theorem 
follows based on (l30l , OTT i. Inequality (fTTT i. and the definition 
of s.c. service curve. ■ 

We have the following theorem for service guarantee. 

Theorem 3: Assume that a data flow A(t) with the traffic 
arrival curve A ^ v b< f, a > is input into a node. Assume 
that the node has an energy-oblivious service curve S ^ sc < 
g,j3 >. Assume that the node is powered with energy supply 
following C(t) ~ sec < fi,ai,f 2 ,a 2 >■ 

• For any time t > 0, x > 0, 

Prob{D{t) > h(a(t) + x, (3 A d 2 (t))} < f ® g(x) (32) 



where Sc is defined by (refeq:SC). This is because if S(t) < 
Sc(t), the energy constraint does not play a role, and the 
actual service amount is equal to S(t). Otherwise, the actual 
service is equal to the amount allowed by the energy con- 
straint, i.e., S c (t). In other words, S e (t) = S(t) A S c (t). 

Next, we prove that S e (t) ~ sc < g,(3 A d 2 >. Denote the 
arrival flow as A(t) and the output flow as A*(t), which is 
equal to S e (t). For any time < s < t, based on Lemma [2] 
and Lemma [3] we have 

A®(/3/\ d 2 )(s)-A*(s) < (A®p(s))Ad 2 (s)-A*(s) 
= (A ® (3(s)) A d 2 (s) - S(s) A S c (s) 
<[A® /3(s) - A*(s)}+ + [S c (s) - d 2 (s)} + 

(28) 



where g is defined with d25l ). d 2 is defined with 
and h(a(t) + x, f3 A d 2 (t)) is the maximum horizontal 
distance between functions a(t) + x and j3 A d 2 (t). 
• For any t > and x > 0, the backlog B(t) is bounded 
by 

Prob{B(t) > x} < f (g) g(x - a (/3 A d 2 )(0)). (33) 

Proof: Theorem|3]follows directly by applying Theorem[2] 
to existing results of service guarantee in traditional stochastic 
network calculus (refer to Chapter 5 of iflOll ). ■ 

C. The Danger of Energy-Oblivious Service 

An argument for using deterministic energy management is 
that the system would be safe if the average energy depletion 
rate equals the predicted energy replenishment rate. Since no 
absolute guarantee can be made on the accuracy of energy 
prediction, there is a chance that an energy outage occurs, i.e., 
the remaining energy-level is lower than the safety threshold. 
In this section, we show that the energy outage probability is 
non-negligible. 

Theorem 4: Assume that a node provides an energy- 
oblivious service curve S ~ S sc< <?,/? >■ Assume that the 
node is powered with energy supply following C(t) ^ sec < 
fl,ai,f 2 , a 2 >. Denote C£(f) = / Q * V(S'(x))dx, and E{t) = 
C(t) — Cg(t), For any time t > 0, given the energy safety 
threshold x > 0, 

Prob{E(t) < x} > 1 - f 2 ® g(x - a 2 /3(0)), (34) 
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where $(t) = J * V{j3{x))dx. 

Proof: For any t > s > 0, we have 

E(t) =C(t)-C* s (t) 

=C{t) - a 2 {t) + 0(f) - C* s (t) + a 2 (t) - /3(f) 
< sup {C(s, t) — a 2 (t — s)} 

0<s<t 

+ sup {P(s,t)-C* s {s,t)} + su V {a 2 (t)-P(t)} 

0<s<t t>0 

= sup {C(s,t) - a 2 (t — s)} 

0<s<t 

+ sup {J3(s, t) - C* s {s, t)} + a 2 $(0). (35) 

0<s<t 

Based on Inequality (JTTJ, the definition of s.e.c. curve, and the 
definition of s.s.c. curve (we note again that the calculations 
of $ and Cg do not change the bounding function), we obtain 
for any x > 0, 

Prob{E(t) > x} < f 2 ® g(x - a 2 /3(0)). (36) 

The theorem follows. ■ 
As an example, assume that there exists an upper curve 
on the energy charging rate of a 2 (t) = pt + a, with a 
bounding function f 2 (x) = e~( x+2 \ where a is a given 
safety threshold energy level and p is the average energy 
charging rate. Assume that the node could provide energy- 
oblivious service S ^ s . c .< 0, g >, where g(x) — e~( x+2 \ 
Assume that 0(t) — J Q P(0(x))dx = pt. Clearly, the average 
energy charging rate is equal to the average energy discharging 
rate and initially the energy level at the node is safe. With 
Theorem H] however, it is easy to calculate that the energy 
outage probability is no less than 72.9%, regardless of the 
time and the values of p and a. 

Remark 8: Theorem [4] and the above example illustrate 
the danger of energy-oblivious service. Since a service is 
practically feasible only if it does not cause an energy outage, 
any service scheduling algorithm should consider the energy 
constraint. This is the main reason why in Section [ill] the 
energy discharging model is inherently coupled with the 
energy charging model and why in our analysis for service 
guarantee, we enforce an energy constraint on the energy- 
oblivious service. 

V. Stochastic Service Guarantee for a 
Networked System 

In the previous section, we analyzed the single node case. 
Here, we study the service guarantee question in a network 
for which the network nodes are powered with a renewable 
energy source. It is well-known that we cannot simply add 
the delay bounds for each individual nodes to obtain the end- 
to-end delay bound along a path 0, lfl2l . since otherwise the 
network-wide delay bounds would be too loose to be useful. 
As a solution, the concatenation property of service curves 
should be proved and applied iflOl . 

Based on Theorem[2]and the concatenation property of s.c. 
service model, the following theorem holds immediately. 

Theorem 5: Consider a traffic flow passing through of a 
network of N nodes in tandem. Assume that each node 
i(= 1,2, . . . , N) would provide an energy-oblivious service 



curve S l ^ sc < g\p l > to its input. Assume that each node 
i is powered with an energy supply following C % ^ sec < 
fl,a\, f 2,a>2 >■ The network guarantees to the flow a stochas- 
tic service curve S ~ sc < <?,/?> with 

.</ il //'••• // A - (37) 

P = P 1 ® /3 2 P N , (38) 

where 

g\x)=g l ®P 2 (x) : (39) 

/3*(i) = /3'(t) A / p-VIO))*^. (40) 
J o 

Remark 9: Theorem [5] indicates that we can treat the con- 
catenation of multiple nodes as a single system. To obtain the 
end-to-end performance, we just need to apply Theorem|5]into 
Theorem [3] that is, replacing g and A d 2 (t) in Theorem [3] 
with g and (3 in Theorem [5] respectively. 

Remark 10: In Sections |IV] and [V] we only used the upper 
curve of the charged energy amount to constrain a service. 
Actually, another set of similar analysis could be done, if we 
use the lower curve of the charged energy amount to constrain 
a service. The practical meaning of this constraint is to avoid 
overcharging the system and to effectively use the energy, i.e., 
the service should deplete the energy at a rate no less than the 
minimal charging rate. 

VI. Nodes Powered with Multiple Energy Sources 

It is possibly that a network node might be charged with 
multiple energy sources. For example, several types of energy 
such as solar, wind, vibrational, and thermal among others 
can be scavenged from the surroundings of a sensor node to 
replenish its battery lfl9l . In this case, we need to study the 
statistical features of the integrated energy source, for which 
we have the following theorem. 

Theorem 6: Assume that a node has N energy sources, 
denoted by C 1 , C 2 , . . . , C N , with & following &{€) - sec < 
f\,ot\, f 2 ,a 2 >. Assume that the hardware permits the multi- 
ple energy sources to charge the battery simultaneously. The 
overall energy supply to the node, denoted by C, follows 
C{t) ~ sec < fi,cei,f 2 ,a 2 >, where 

fi(x) = fl®fl ■ ■ ■ ®f?{x) -(N- 1), (41) 
ai{t) =a{+al... + a?(t), (42) 

h(?) = /2 «"/!•■• ®f?(x), (43) 

a2 {t)=a\+al... + a%{t), (44) 

Proof: We only prove Equations (HTl i and (l42l by apply- 
ing Inequality (1121 1. as the proof for Equations d43l and (l44l 
is similar by using Inequality (fTTl i. In addition, we only need 
to prove the case that N = 2, because the theorem holds by 
recursively applying the result for N = 2. 



As we assume that the multiple energy sources can charge 
the node simultaneously, for any time t > 0, C(t) — C 1 (t) + 
C 2 (t). For any t > s > 0, we have 

inf {C(s,t) -ai(i-s)} 

0<s<t 

= inf {C 1 (s,t) + C 2 (s,t)-(a 1 1 (t-s) + a 2 1 (t-s))} 

0<s<t 

= inf {C 1 (s,t)-a 1 1 {t-s) + (C 2 (s,t)-a 2 1 (t-s))} 

0<s<t 

> ini : {C\s,t) -a\(t-s)+ ini : \C 2 (s,t) - a\{t - s)]} 
= ini I {C\s,t) -a\(t- s)} + M JC 2 ( S ,t) - a 2 (t - s)} 

0<s<t 0<s<t 



(45) 

The lower curve and the lower bounding function of C(t), i.e., 
Equations ( |4TT > and (|42] |. are thus proved based on the above 
inequality and Inequality dl2l . ■ 

VII. Further Improvement on Bounds 

So far, we have developed an analytical framework based 
on stochastic network calculus to evaluate the performance of 
network systems with nodes powered by renewable energy 
sources. This analytical framework is general enough no 
matter whether or not the energy-oblivious service process and 
the energy charging process are independent, and no matter 
whether or not multiple energy sources, if they exists, are 
independent. If we assume that the energy-oblivious service 
process and the energy charging process are independent or 
that multiple energy sources are independent, which we believe 
is true in most applications, better performance bounds could 
be obtained, using the following lemma. 

Lemma 5: Assume that the complementary cumulative dis- 
tribution functions ( CCDF) of non-negative random variables 
X, Y are Fx £ T and Fy e T, respectively. Assume that 
fi(x) < F x {x)_ < f 2 (x), 5i (x) < F Y {x) < g 2 (x), and 
/ij/2,51,32 G J~. Denote Z = X + Y. Assume that X and 
Y are independent, there holds for V.x > 0, 



and 



Fz(x) <l-/2*32(x), 



Fz(x) >l-A*gi(x), 



(46) 



(47) 



where / x = 1 - [/i]i,/ 2 = 1 - [/ 2 ]i,0i = 1 - Mi, 02 = 
1 — [02] i, and * is the Stieltjes convolution operation [17]. 

With Lemma [5] if the independence assumption could be 
made, the bounds in the theorems of this paper could be 
revised accordingly. That is, we could apply Lemma [5] instead 
of Lemma Q] in the proofs of the theorems to get better 
bounds. For example, if we assume that two energy sources 
are independent, the lower and upper bounds in Theorem [6] 
should be revised, respectively, to 



/ 1 (x) = l-(l-/ 1 1 )*(l-/ 1 2 ), 
/aOc) = !-(!- /!)*(l-/f). 



(48) 
(49) 



As an example, assume that fl(x) — f 2 {x) = e~ 2x and 
f£( x ) = fi( x ) — e ~ x - With Theorem [6] we have the lower 
bound fi(x) = e~ 2x and the upper bound f2(x) — 2e~ x ^ 2 . 
But, by applying Lemma [5] we can get new lower bound 
fi{x) = (1 + 2x)e~ 2x and new upper bound f2(x) = 




Fig. 1. Comparison of old and new bounds 

(1 + x)e~ x . It is easy to verify that the new bounds are tighter 
as shown in Fig. [TJ 

VIII. Related Work 

With the increasing importance of greening computing, QoS 
and performance evaluation of network systems with uncertain 
energy supply have attracted much attention. Many perfor- 
mance modeling approaches have been developed, which 
could be divided into two main categories, deterministic and 
stochastic approaches. 

In the deterministic group, Zafer l23l et al. use deterministic 
network calculus to model traffic arrival and traffic departure. 
They use a power-rate function to link the traffic departure 
rate and energy consumption rate. By considering the special 
features of specific power-rate functions, they formulate and 
solve the optimal transmission scheduling problem under the 
given energy constraints. Their work only focuses on single- 
node analysis and assumes that traffic arrivals and service 
rates are deterministic. Kansal et al. ifTTl propose a so-termed 
harvesting theory to help energy management of a sensor 
node and determine the performance levels that the sensor 
node can support. The basic idea of the harvesting theory 
is to use a leaky bucket model to represent energy supply 
and energy depletion. Moser et al. lfT31l describe energy-aware 
scheduling and prove the conditions for a scheduling algorithm 
to be optimal in a system whose energy storage is replenished 
predictably. 

In the stochastic group, Markov chain models have been 
used extensively. Susu et al. [20 1 use a discrete-time Markov 
chain in which states represent different energy levels. Some 
work [16 1 uses a Markov chain model to capture the influence 
of clouds and wind on solar radiation intensity. Relevant to 
stochastic energy modeling, there are many efforts to predict 
a stochastic energy supply. Lu et al. lfT4l assess three pre- 
diction techniques: regression analysis, moving average, and 
exponential smoothing. Recas et al. lfl8l propose a weather- 
conditioned moving average (WCMA) model, which adapts to 
long-term seasonal changes and short-term sudden weather 
changes. Moser et al. lfT31l introduce energy variability curves 
to predict the power provided by a harvesting unit. 

We develop our analytical framework based on stochastic 
network calculus 0, 0, J9j , iflOll . Unlike deterministic net- 
work calculus JJJ, lfl2l . which searches for the worst-case per- 
formance bounds, stochastic network calculus tries to derive 
tighter performance bounds, but with a small probability that 
the bounds may not hold true. Since most renewable energy 



s 



sources, such as solar and wind energy, are not deterministic, 
stochastic network calculus is a good fit for the performance 
evaluation of systems using renewable energy. Nevertheless, 
traditional stochastic network calculus was not originally 
targeted at modeling such systems. Substantial work is thus 
required to extend this useful theory. 

Recent interesting work by Wang et al. [22 1 uses stochastic 
network calculus to evaluate the reliability of the power grid 
with respect to renewable energy. Their energy supply and 
demand models are a subset of the models we present in 
Section Hn] Their work shows a good example of how to tailor 
our models for a specific application. Another fundamental 
difference is that Wang et al. define the energy supply and 
energy demand as two de-coupled random processes, while in 
our work energy discharging is inherently coupled with energy 
charging. 

Finally, related to analytical frameworks for performance 
modeling, there is a large body of research on energy-aware 
scheduling algorithms. For example, Niyato et al. lfl6ll in- 
vestigate the impact of different sleep and wake-up strategies 
on data communication among solar-powered wireless nodes. 
In ||2T1| . Vigorito et al. propose an adaptive duty-cycling 
algorithm that ensures operational power levels at wireless 
sensor nodes regardless of changing environmental conditions. 
In Q, Gorlatova et al. measure the energy availability in 
indoor environment and based on the measurement results 
they develop algorithms to determine energy allocation in 
systems with predictable energy inputs and in systems where 
energy inputs are stochastic. In the stochastic model, they 
assume that energy inputs are i.i.d. random variables. Unlike 
the above work, our analytical framework is generic and 
uses only abstract notations on energy charging/discharging 
amount, traffic arrival amount, etc. As such, all the above 
work could be treated as a special case of our more general 
framework in which the abstract functions are replaced with 
concrete ones for specific applications. 

IX. Conclusion and Future Work 

The application of renewable energy poses many chal- 
lenging questions when it comes to a QoS guarantee for 
IT infrastructure. Particularly, the strong call for a generic 
analytical tool for the performance modeling and evaluation of 
such systems remains unanswered. We, in this paper, develop 
such a theory using a stochastic network calculus 0, |6|, 
1 9 1, 1 10 1 approach. We enrich the existing stochastic network 
calculus theory to make it useful for evaluating systems with 
a stochastic energy supply. We introduce new models to 
capture the dynamics in the energy charging and discharging 
processes, and derive new analytical results to provide funda- 
mental performance bounds, which could be used to guide 
the energy management and the design of task scheduling 
algorithms. 

This paper mainly focuses on the introduction of a theo- 
retical framework. Along the line, we envisage many future 
research topics, including (1) the application of energy charg- 
ing/discharging models in different applications (e.g., |20|, 
1 22 1), (2) analyzing the service model of different scheduling 
strategies (e.g., [15]), and (3) network reliability analysis (e.g., 
network-wide energy outage analysis). 
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Appendix 

A. Proof of Lemma [7] 

Proof: The proof of (fTTT i could be found at |9 j. We only 
prove dT2b - For any y > x > 0, 



{X + Y < y} n {X > x} n {Y > y - x} = 0, 
where denotes the null set. We thus have 

{X + Y > y} D ({X > x} n {Y > y - x}) 
and hence 

Prob{X + Y > y}> Prob{{X > x} Ci {Y > y ■ 
>1 - Prob{X <x}- Prob{Y < y - x} 
=Prob{X > x} + Prob{Y > y - x} - 1. 



(50) 
(51) 

?}} 

(52) 



Since the above inequality holds for all y > x > 0, we get 

Prob{X + Y >y}> sup {Prob{X > x] 

y>x>0 

Prob{Y > y - x}} - 1, (53) 



Case 1: /i(s x ) + f 2 (t - 8l ) < f 3 (t). The right hand 
side of Inequality (IBTt equals /i(si) + / 2 (t — si). Since 
he J 7 and /i(0) = 0, we have / x (si) > 0. Therefore, 

Ht- Sl )<f 3 (t). (59) 

We need to show that the left hand side of Inequality ( |57| > 
is no larger than /i(si) si). For the given t > 0, 

there exists a S2(0 < S2 < t) such that 

inf {/i(«)+/ 2 (t-s)A/ 3 (t)} =/ 1 (« 2 )+/ 2 (*-*2)A/ 3 (*). 

0<s<i 

(60) 

We have 

/i(s 2 ) + h(t - s 2 ) A f 3 (t) < h(si) + h(t - .si) A h{t) 

=h(si) + h(t- Sl ). 

(61) 

In the above, the first inequality is due to (1601 . and the 
second equality is due to ( |59l . 

Case 2: /i(si) + /2(i — si) > h(t). The right hand side 
of Inequality d57| > equals h(t). We also have /i(0) + 
h(t - 0) > /i(si) + / 2 (t - si) due to d58]l. Therefore, 



/2(t)>/i(*i) + /2(t-ai)>/ 3 (*), 



and 



which is 



Fz(y) > Fx®F Y {y) - 1 > h®9i{y) - 1- 



(54) 



inf {h{s) + h{t -a) A hit)} < h(0) + hit) A hit) 

0<s<t 

= fs(t) 



■ (62) 

B. Proof of Lemma\2\ In conclusion, Inequality (ED is true for both cases. The 

Proof: By definitions of the ® and A operations, we need lemma is proved. ■ 

to prove that „ „ „ r , rgi 

C. Proof of Lemma \3\ 

inf {/i(s)+/ 2 A/ 3 (t-s)} < inf {/i(s)+/ 2 (t-s)}A/ 3 (t). f ro0 /:- To evaluate the value of X x A X A - X 2 A X 3 , we 

(55) have: 

Since / 3 e J 7 , we have 



inf {his)+hAf 3 (ts)} < inf {/i(«)+/ 2 (t-s)A/ 3 (t)}. 

0<s<i 0<s<t 

(56) 

We therefore only need to prove that 

inf {A( S )+/ 2 (t- S )A/ 3 (t)} < ini M/ 1 ( s )+/ 2 (t_ s )}A/ 3 (i). 

0<s<t 0<s<t 

(57) 

First, for any given t > 0, there exists a si(0 < si < t), such 
that 

inf {his) + hit - s)} = hist) + hit - si). (58) 

0<s<t 

In addition, when t is given, hit) is a constant. 

We next prove that Inequality (|57| i holds true for the 
following two exclusive cases: 



Case 1: X t < X A and X 2 > X 3 . We have = X 2 = 
X 3 = X 4 , because X x > X 2 > and X 3 > X 4 > 0. 
Thus (EHi holds. 

Case 2: Xi < X A and X 2 < X 3 . We have 

X\ A X4 — X 2 A X 3 = Xi — X 2 < X\ — X 2 + X 3 — X A . 

Case 3: Xi> Xi and A" 2 > X 3 . We have 

Xi A X4 — X 2 A X 3 — X4 — X 3 < 0. 

Inequality d24l > holds since the right side is no less than 
0. 

. Case 4: X x > X 4 and X 2 < X 3 . We have 

X\ A X4 — X 2 A X 3 ^ X\ — X 2 < ^1 — A^ 2 — J^ 3 — X4 . 

The lemma is proved since the above list covers all possible 
scenarios. ■ 



10 



D. Proof of Lemma^} and F Y < gi- Hence, we have 

Proof: The proof of d46l could be found at fT71 . We only f x 

prove <g7j. ' Fz(x) = J F x (x - y)dF Y (y) < j ^ fi(x - y)dF Y (y) 

For independent non-negative random variables, X and Y, „x 

we have Va; > 0, = / F Y (x — y)dfx(y) 



F z (x) = F x (x - y)dF Y (y). < / _ y)d/i(y) = f^- gi{x) _ (63) 



Note that F x ,F Y Ji, gi are wide-sense increasing, F x < A Note mat the third equa i ity ho lds because Stieltjes convolution 

operation is commutative. Inequality ( f4Tb is thus proved since 
F z (x) = l-F z (x). ' M 



